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In this paper we investigate the thermal evolution of hybrid stars, objects composed of a quark 
matter core, enveloped by ordinary hadronic matter. Our purpose is to investigate how important 
are the microscopic properties of the quark core to the thermal evolution of the star. In order to 
do that we use a simple MIT bag model for the quark core, and a relativistic mean field model for 
the hadronic envelope. By choosing different values for the microscopic parameters (bag constant, 
strange quark mass, strong coupling constant) we obtain hybrid stars with different quark core 
properties. We also consider the possibility of color superconductivity in the quark core. With this 
simple approach, we have found a set of microscopic parameters that lead to a good agreement with 
observed cooling neutron stars. Our results can be used to obtain clues regarding the properties of 
the quark core in hybrid stars, and can be used to refine more sophisticated models for the equation 
of state of quark matter. 
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I. INTRODUCTION 



Scientists have been trying to determine the equation 
of state and composition of compact stars for a long time 
([U, and references therein). Many of these studies 
determine constraints for the equation of state of stel- 
lar matter so that the predicted stellar properties agree 
with the observed masses and radii of pulsars. One of the 
biggest practical shortcomings of such techniques is the 
lack of accurate and reliable observations for the radii of 
compact stars, which makes constrainingthe equation of 
state very challenging. A recent paper Q has put forth 
unprecedented constraints on the mass-radius diagrams 
of pulsars, which would certainly make the range of pos- 
sible equations of state narrower. The analysis of this 
data however has been questioned and therefore is 
not universally accepted. Another method that might be 
used to complement the results provided by the analysis 
of the structure of the compact star is the examination 
of the thermal evolution of the object. The cooling of a 
compact star is intrinsically connected to its microscopic 
composition, and therefore different microscopic models 
might lead to different thermal evolutions. The results 
of the theoretical investigations might then be compared 
with the wealth of observed data of the thermal proper- 
ties of compact stars @, Q . In this paper we will follow 
the thermal evolution approach, and with this probe the 
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composition of compact stars. 

In the following we will consider the thermal evolu- 
tion of hybrid stars (HS), which are objects composed 
of a strange quark matter core, enveloped by ordinary 
hadronic matter [U, In this work, neutron stars (NS) 
denote objects made up solely of hadronic matter. 

The cooling of neutron stars is dominated by neutrino 
emission for the first 1000 years (possibly more for the 
slow cooling scenario) , being replaced by surface pho- 
ton emission after that. Amongst the neutrino emission 
mechanisms, the direct Urea (DU hereafter) process sets 
itself apart by being one of the most efficient ones, with 
emissivity of the order of 10^^ erg cm~^ s~^ Due to 
momentum conservation, the DU process can only take 
place if the proton fraction is above a certain value. The 
actual proton-fraction, above which the DU process is 
allowed to take place, depends on the underlying micro- 
scopic model, and it usually lies between 11%— 15%. If 
a neutron star allows the DU process to take place, the 
star will exhibit a fast cooling (sometimes also referred 
to as enhanced cooling). If neutron superfluidity [Tol.[Tl| 
(superconductivity in the case of protons) occurs, pairing 
will suppress the emission of neutrinos from hadronic pro- 
cesses and smooth the so-called enhanced cooling. The 
situation for hybrid stars is significantly different, since 
the presence of deconfined quark matter plays an impor- 
tant role for the cooling of the object. The coolin g of hy - 
brid stars has been addressed by many authors lol .jlSl 14 1 . 
In this work, differently than in references 0, [ij, fl4 |. 
where the cooling of hybrid stars with the same equation 
of state was studied, we will examine the cooling of hy- 
brid stars with different EoS's for the quark phase. As a 
result, we will obtain quark cores with different proper- 
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ties (mass, radius, composition). This should provide us 
with a measure of how important such properties are for 
the thermal evolution of hybrid stars, which might aid us 
in constraining the EoS of compact stars. In order to per- 
form such study we will use a MIT bag model equation 
of state for the quark phase, and a relativistic-mean-field 
model EoS for the hadronic phase. The different EoS's for 
the quark phase arc obtained by varying the microscopic 
parameters of quark matter, namely the bag constant B, 
the strong interaction constant , and the strange quark 
mass TTig. The macroscopic structure of the star is ob- 
tained by solving the Tolman-Oppenheimer-Volkoff equa- 
tions [13, [13 ■ We note that there are more so phi sticated 
models for hybrid stars, see for example [19l - l2ll |. How- 
ever, the relatively simple MIT bag model is appropriate 
for our study. Due to its phenomenological nature, it per- 
mits one to change its microscopic properties with ease, 
allowing us to obtain quark cores with different proper- 
ties. This in turn allows us to investigate how such quark 
cores effectively influence the thermal evolution of hybrid 
stars. The aim of this procedure is to obtain constraints 
on the properties of cold quark matter, which then can 
be used in more sophisticated microscopic calculations. 

This paper has the following structure: in Sec. |TT] we 
describe the microscopic model used, and present the re- 
sults for the macroscopic structure of the stars; in Sec. lIIII 
our results for the thermal evolution of the stars discussed 
in Sec. In] are presented; in Sec. IIVI the role of supercon- 
ductivity in the quark core is discussed, and the results 
for the cooling of superconducting hybrid stars are pre- 
sented; our final remarks and conclusions can be found 
in Sec. |V| 



II. MODEL 

We use a standard relativistic mean field model [l|,[H to 
describe the hadronic phase, in which the interactions be- 
tween the baryons are mediated by meson fields (cr, w, p). 
The parametrization used is G300, which can success fully 
reproduce the properties of saturated nuclear matter [23] ■ 
We have also allowed pairing in the hadronic phase. The 
pairing patterns included are the neutron singlet (^S'o) 
and triplet {^P2) states, as described in [l^, It is 
also possible that protons can form pairs, becoming a 
singlet superconductor; however, the pairing of protons 
in the core is still not very well understood (see [12l | and 
references therein). For this reason have decided to take 
a conservative approach and only consider neutron su- 
perfluidity. 

For the quark matter description we use the MIT Bag 
Model, with first order corrections in the strong interac- 
tion coupling constant (as) [23l - l26| . In the framework of 
the MIT Bag Model, the quarks present at the relevant 
densities are the up, down and strange quarks. It is im- 
portant to note that, since we are considering massive 
strange quarks, electrons are required for charge neutral- 
ity to be achieved. At lower densities, the electron popu- 



TABLE I: Different parameter sets used to obtain the quark 
matter EoS. 



Parameter Set 


fii/^ (MeV) 


m, (MeV) 


as 


A 


139 


200 


0.7 


B 


150 


175 


0.4 


C 


160 


150 


0.1 



lation is more significant, due to suppression of the mas- 
sive strange quarks. This is an important phenomenon, 
that may lead to the formation of ultra-stron g el ectric 
fields near the surface of bare strange stars j27H3g | . 

In this work we also take into account the possibility 
of color superconductivity. The pattern that will be con- 
sidered is the Color-Flavor-Locked phase [3l|, where all 
quarks are paired. The CFL phase is the most likely con- 
densation pattern at densities > 2po (32j . For intermedi- 
ate densities (^ 2/9o) model calculations seem to indicate 
that quark matter is in a 2SC phase [32|. Another possi- 
bility is that quark matter forms a crystalline supercon- 
ductor, where the momenta of the quark pairs do not add 
to zero j33l [s^ l . Since we are considering quark matter 
only at the inner core of compact stars, we will consider 
only the CFL phase. It should be noted that one expects 
corrections to the quark matter EoS if pairing is consid- 
ered. As shown in tSa, such corrections are ^ A^/i^. 
The authors in ref. [361 ] have shown, however, that the 
effects of such corrections to the structure of the star are 
only noticeable for A > 50 MeV. Therefore, for the val- 
ues of A considered in this study (0.1 — 10 MeV) they 
can be safely disregarded. Note that as the small values 
for the gap (A) do not have any appreciable effect on the 
equation of state, our analysis is still valid for any quark 
pairing scheme (not necessarily color superconductivity) , 
as long as it affects all quark flavors in similar way. In 
addition to the quark matter and hadronic matter equa- 
tion of state, we have also included a crustal EoS, which 
is given by the Baym-Pethick-Sutherland EoS [37| . 

The equation of state for the compact star is calculated 
assuming local charge neutrality and beta equilibrium, 
as is usual for any bulk system. The phase transition 
to quark matter is also calculated in the usual way, by 
connecting the two EoS's at the point where the quark 
EoS has a higher pressure. The different values used 
for the quark EoS parameters are illustrated in Table [J 
These parameters were chosen to cover a wide range of 
quark core properties (mass, radius and composition), 
which might potentially lead to different thermal evolu- 
tion. The final EoS's obtained are shown in Fig.[Tl where 
we plot the pressure as a function of energy density for 
the different EoS's investigated. 

The mass-radius diagram for each EoS studied can be 
found in Fig. [2l The onset of quark matter can be easily 
identified as the sharp kink found in the sequences. Note 
that for simplicity, we adopt a Maxwell construction for 
the phase transition. We do not expect any qualitative 
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FIG. 1: (Color online) Equation of state for the hadronic and 
quark phases shown for different parameter sets. 
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FIG. 3: (Color online) Mass-radius diagram for the quark core 
of the hybrid stars shown in Fig. (2] Also shown the quark core 
radius for the stars used in the cooling simulations (see table 
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III. RESULTS 
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FIG. 2: (Color online) Mass-radius diagram for stars whose 
equations of state are shown in Fig. [T] The sharp kinks denote 
the phase transition to quark matter. 



change of our conclusions if a more elaborate construction 
were used. 

Given that each of the sequences shown in Fig. [2] has 
a different quark matter EoS, the quark core of the stars 
in each sequence should exhibit different properties. This 
is displayed in Fig. [3] where the mass-radius diagram for 
the inner quark core of the hybrid stars is shown. 

Fig. |3] clearly shows the distinct properties of the quark 
core for the different EoS's used. We now turn our at- 
tention to the thermal evolution of the stars shown in 
Figs. [5] and [31 and determine the role that the different 
microscopic compositions play in the thermal evolution 
of the object. 



The cooling of compact stars is given by the thermal 
balance and thermal energy transport equation (G = c = 



5(Ze^ 
dm 

djTe'f') 
dm 



P\/\ — 2m/r 



e.e^^ + c,,^^") ,(1) 



dt 



Wir'^r^Kpy^l — 2m/ r 



(2) 



In Eqs. |T1 and |2l the structure of the star is given by 
the variables r, p(r) and m{r), that represent the radial 
distance, the energy density, and the stellar mass, re- 
spectively. The thermal variables are given by the tem- 
perature T(r,t), luminosity l{r,t), neutrino emissivity 
e^{r,T), thermal conductivity K{r,T) and specific heat 
Cy{r,T). The boundary conditions of ([T]) and ^ are 
determined by the luminosity at the stellar center and 
at the surface. The luminosity vanishes at the stellar 
center since at this point the heat flux is zero by defini- 
tion. At the surface, the luminosity is defined by the re- 
lationship between the mantle temperature and the tem- 
perature outside of the star 0, S IHj 113 ■ The micro- 
scopic inputs that enter Eqs. ^ and ^ are the neutrino 
emissivities, specific heat and thermal conductivity. For 
the hadronic phase, the neutrino emission processes con- 
sidered in this work are the direct Urea, modified Urea 
and bremsstrahlung processes. A detailed review of the 
emissivities of such processes can be found in reference 
[40j . We note that the neutrino emission involving neu- 
trons, and the neutron specific heat are exponentially 
suppressed as the temperature drops below the neutron 
superfluidity critical temperature. A description of the 
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TABLE II: Stellar properties of the stars used for the thermal 
evolution investigation. All symbols have their usual mean- 
ing, and Mq and Rq denotes the mass and radius of the quark 
matter core, respectively. 



Parameter set 


Mq/M 


-R, (km) 


R (km) 


M (Mq) 


A 


0.48 


8.28 


13.26 


1.33 




0.60 


9.14 


13.15 


1.55 


B 


0.62 


8.21 


12.20 


1.33 




0.74 


8.86 


11.82 


1.55 


C 


0.82 


8.21 


10.89 


1.33 




0.89 


8.29 


10.03 


1.55 



shown Figs. [4] and [5l The reason for such behavior is the 
presence of the hadronic superfluidity, which suppresses 
the hadronic DU process. Without superfluidity a neu- 
tron star with higher mass could exhibit a faster cooling 
than a lower mass star, since the DU process could be 
active in the former but not in the latter. The fact that 
the hadronic layer plays the same role, regardless of the 
total mass of the object, is very important for our anal- 
ysis, since it allows us to better understand the effects 
introduced by the different quark cores in the thermal 

p-\rnhitinn nf thp star 



paring model, and critical temperature can be found in 
reference As for the quark core, we consider the 

equivalent processes involving quarks. These are the 
quark direct Urea (QDU), quark modified Urea (QMU), 
and quark bremsstrahlung processes (QBM). The emis- 
sivities of such processes were calculated in [4l| . The spe- 
cific heat of the hadrons in the hadronic phase is given 
by the usual specific heat of fermions, as described in 
For the quark phase we use the specific heat as calcu- 
lated in [4l[ . Finally for the thermal conductivity of the 
hadronic matter, we follow the calculations of (42|, and 
for the quark matter we use the results of fisj . 

In order to investigate how the difference in the quark 
cores manifests itself in the thermal evolution, we inves- 
tigate the cooling of hybrid stars of same mass, but dif- 
ferent quark cores. The properties of such stars are listed 
in Table [ni 

The thermal evolution of the objects from table |lT] is 
illustrated in Figs. [4] and El which shows the redshifted 
temperature as a function of the age of the star. We have 
also plotted prominent observed temperatures for a few 
pulsars. The data is separated into two sets. The first 
set, labeled spin-down, contains the temperature of ob- 
jects whose age are estimated based on their spin-down 
age. The second set, labeled kinematic, consists of age 
estimates based on the observed motion of the object 
with respect to its progenitor supernova remnant. A de- 
tailed discussion about the observed data can be found 
in references d, @ . 

The results indicate that although there is little differ- 
ence between the cooling of stars with different gravita- 
tional masses, we can see a significant difference between 
stars with different quark core properties. The differ- 
ent cooling tracks exhibited by stars of different models 
(A. B and C) can be attributed to the differences in the 
quark core. Figs. 3] and [5] show that for a more mas- 
sive quark core, the cooling of the star becomes faster. 
This is a result of the microscopic composition of such 
models, that allow for more intense presence of the QDU 
(due to a higher electron population). It is also impor- 
tant to understand that regardless of the total mass of 
the star, the hadronic layer plays the same role. This 
is indicated by the similar cooling behavior shown for 
the purely hadronic stars of masses 1.33 and 1.55 M© 




Age (years) 

FIG. 4: (Color online) Cooling of hybrid stars with gravita- 
tional mass of 1.33 Mq. Too denotes the temperature observed 
at infinity, and the x-axis the age in years. Also plotted is the 
cooling of the equivalent purely hadronic star with the same 
mass. Two sets of observed data were used, one whose age es- 
timate is based on spin-down rate, and another that is based 
on the kinematic age [1, Q . 



IV. ROLE OF COLOR SUPERCONDUCTIVITY 

It seems clear that the presence of the quark core en- 
hances the cooling of the star to a point that is it in 
clear disagreement with the observed data, whereas the 
purely hadronic model seems to do a (slightly) better 
job. We note however that results presented so far do 
not represent the whole picture. As mentioned in section 
im strange quark matter is expected to be in a super- 
conducting phase. The most likely condensation pattern 
for high densities is the Color-Flavor-Locked phase, in 
which all quarks are paired. Because of pairing the di- 
rect Urea process is suppressed by a factor e~'^^'^ , and 
the modified Urea and the Bremsstrahlung process by 
a factor e"^'^/"^, for T < Tc, where A is the gap pa- 
rameter for the CFL phase and Tc is the critical tem- 
perature below which strange matter undergoes a phase 
transition into CFL matter. In addition to changes 
in the neutrino cmissivities, the specific heat of quark 
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FIG. 5: (Color online) Same as Figure 14] but for stellar mass 
of 1.55 Mq. 
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FIG. 6: (Color online) Cooling of hybrid stars with gravita- 
tional mass of 1.55 Mq, and A = 0.1 MeV. Axis and observed 
data are the same as in Fig. 2] 



matter is also modified, being reduced by the factor 
3.2(rjT)*(2.5-1.7(T/rc) + 3.6(T/Tj2)e-^/T [|. The 

critical temperature for the CFL phase is currently not 
known, although it is believed to be smaller than the 
standard Bardeen-Cooper-SchriefFer {Tc ~ 0.57A), due 
to instanton— anti-instanton effects. In this work we fol- 
low the footsteps of the authors in and use Tc = 0.4A. 

In addition to the traditional neutrino cmissivities, 
one could also include the cmissivities from massive and 
massless Goldstonc modes. These emissions, as discussed 
in [l^, are small, and should have small effects on the 
thermal evolution of the object if more powerful emission 
processes are present. In addition to that, we have also 
not considered the contributions of photons, and Gold- 
stone modes to the thermal conductivity of CFL quark 
matter [l^. If these modes are considered, the thermal 
conductivity of quark matter should be higher than that 
calculated in [4J]. Not denying that such contributions 
are important, we do not expect that they would strongly 
affect our results. As discussed in [l6[, a higher thermal 
conductivity would only mean that the quark core will 
be in thermal equilibrium with the hadronic phase more 
quickly. We intend to include these contributions in fu- 
ture calculations. 

We now show in Figs.[S]and[7]the results for the cooling 
of hybrid stars of 1.55 Mq, whose quark core is composed 
of strange quark matter in the CFL phase with A = 0.1 
MeV and A = 1.0 MeV, respectively. 

As discussed in the last section, the hadronic layer 
plays a small role in the cooling of hybrid stars, with the 
quark core dominating the thermal evolution (at least 
for the first 1000 years). Following this logic it is only 
obvious that the introduction of the CFL phase should 
alter the cooling substantially. This is indeed the case, as 
shown in Figs. [S] and [71 We see that objects with A = 0.1 
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FIG. 7: (Color online) Same as Fig. [7] but for A = 1.0 MeV 
MeV. 



MeV exhibit a much slower cooling than their unpaired 
counterparts ( Figs. [4] and [5|). The slower cooling is due 
the suppression of the neutrino emission processes in the 
quark core resulting from quark pairing. Fig. [7| shows 
that a higher value of A will result in even slower cooling 
as a higher A will lead to stronger suppression. We note, 
however, that for A > 1.0 the resulting thermal evolu- 
tion is almost the same as that of A = 1.0. This is due 
to the fact that the exponential exp(— A/T) effectively 
saturates for A > 1.0 MeV. 

Figs. Eland [7] also show that hybrid stars, whose quark 
cores are in a CFL phase, agree relatively well with the 
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served data can be better explained by hybrid stars whose 
quark phase is in a CFL state with pairing gaps between 
0.25 and 0.5 MeV. These results also show us that the dif- 
ference between models B and C is very subtle, with both 
these models describing the same data points. Model A, 
however, is substantially different, as can be seen in Fig.[51 
One can see that as opposed to B and C, in the range 
of age between 10^ to 10^ years, the most massive star 
in Model A exhibits a slower cooling. This behaviour is 
due to less massive quark cores in this model. We note 
that there are still two objects that do not fall within 
the range of the cooling curves, these are the two oldest 
stars in the graph. We believe that this might be due to 
uncertainties in the estimate of their spin-down age. As 
pointed out in the spin-down age may be an over- 
estimation of their actual age. If indeed that is the case 
these objects may be younger than they appear, in which 
case they may fall into the range defined by the cooling 
curves. The analysis of the age estimate of these objects, 
however, is beyond the scope of this paper. 



observed data. We see, however, that smaller values of A 
0.1 MeV) lead to agreement with only a few objects 
(those with lower temperatures), while A ^ 1.0 MeV 
leads to agreement with those with higher temperature. 
This behavior hints that the best agreement will possibly 
be obtained for intermediate values of A, if one means to 
interpret the whole set of observed data as hybrid stars. 
In order to investigate this we explore different values of 
A, and calculate the cooling band spanned by the cooling 
curves of stars with 1.00 Mq and the maximum mass of 
each respective parameter set. This result is shown in 

Figs.mni 

The results shown in Figs. [511101 indicate that the ob- 



V. DISCUSSION & CONCLUSIONS 

It is the purpose of this paper to investigate the im- 
portance of the quark core properties to the cooling of 
hybrid stars. Although there have been a handful of 
papers analyzing different models for hybrid stars, and 
their agreement with observational mass and radius con- 
straints [H, [11 [13, m El 0, we intend to provide a 
complementary study, in which observed temperatures 
of pulsars, together with cooling simulations of hybrid 
stars are used to infer the importance of the microscopic 
properties of the quark core to the thermal evolution. 
To achieve that goal, we have used a MIT bag model to 
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describe the quark phase, and a relativistic mean field 
model for the hadron phase. Different values for micro- 
scopic parameters of the quark phase were used, which 
led to quark cores with different properties. 

The results obtained indicate that it is possible to ex- 
plain the observed thermal properties of pulsars with 
a hybrid star model. Furthermore, the best agreement 
with the observed data was obtained for hybrid stars 
with a CFL quark core with as = 0.4, B^^^ = 150 
MeV, ms = 175 MeV and with A = 0.25 - 1.0 MeV. 
We note that parameter set A (a^ = 0.7, S^/** = 139 
MeV, TUs = 200) was also in good agreement with the 
observed data. Furthermore, this parameter set has the 
advantage of agreeing with the observed mass of J1614- 
2230 (1.96 ± O.O4M0) Parameter set A assumes, 

however, that as = 0.7 which might be considered too 
high for the first order corrections used in the calcula- 
tion of the quark matter EoS. Because of this, the re- 
sults obtained from parameter set A should be regarded 
with care, which is the reason why we selected parame- 
ter set B (which also exhibits satisfactory agreement with 
cooling data, albeit not with the mass of J1614-2230) as 
that with the best agreement with the observed cooling 
data (given its more reliable microscopic model) . We can 
however, learn something from the results of parameter 
set A. i.e. that the cooling will present a better agree- 
ment with the observed data for models that yield less 
massive quark cores, and lower electron population (as in 
parameter set A). For higher values of the gap parameter 
(A > 10 MeV), the band of possible cooling curves be- 
comes too narrow, which worsens the agreement with the 
observed data. We note that this result is independent 
of the equation of state used for the hadronic matter, as 
long as the latter presents a suppression of fast neutrino 
emission processes. This is very reasonable, considering 
that the hadronic layer is constrained to the lower den- 
sity regime, in which neutron superfluidity is very likely 
to take place. 

It is important to remember that in this paper we have 
considered a somewhat simple microscopic model, and a 
more sophisticated study, considering for example differ- 
ent degrees of freedom, and chiral symmetry restoration 
[20I [21I l46j would certainly be of interest. Work in that 
direction is currently in progress. We believe, however, 
that even for a more sophisticated microscopic model our 
conclusions should still hold. If so, this would imply that 
a quark core with the properties similar to those of the 
best case scenario (see parameters above) should provide 
the best agreement with the observed data. 

We note that, although we have considered neutron 
superfluidity, we have not included proton superconduc- 



tivity in the core, nor the emission of neutrinos due 
to the pair breaking/formation process. The aforemen- 
tioned processes were the subject of recent papers [islligj 
aimed at modeling the observed temperature evolution of 
the neutron star in Cas A. The inclusion of proton su- 
perconductivity is still speculative, and since we do not 
have a clear picture of how (or if) it manifests itself in 
a high-density environment (see ref. [l2j and references 
therein), we chose to take a conservative approach and 
only consider neutron superfluidity. One could also in- 
clude the pair breaking/formation process, however, this 
should have little effect on our results. Differently from 
references [3 13 j '^e allow fast neutrino processes to 
take place (hence the dip in the cooling curve at around 
100 years), which (even suppressed by pairing) dominate 
over the other neutrino emission processes. Therefore, 
by allowing fast emission processes in the core, we can- 
not describe the observed behaviour of the neutron star 
in Cas A. We note that the cooling of high mass stars of 
parameter set A (in a CFL phase) present temperatures 
and ages comparable to those of Cas A, but the slope of 
the temperature evolution is wrong. This indicates that 
either the neutron star in Cas A has a mass smaller than 
the direct Urea threshold (and thus has no fast neutrino 
emission processes), or that our model needs further re- 
finements. We should also point out that the authors of 
[sot have shown that a hybrid star model can be used 
to model the observed behavior of Cas A, indicating that 
the recent observations do not rule out hybrid stars. Fur- 
thermore, it has been recently proposed that the behavior 
of Cas A can be explained as the late onset of the Direct 
Urea process in spinning down neutron stars (sij . We 
note that the same principle could be applied to hybrid 
stars, where the transition to quark matter could be de- 
layed due to spin-down. This study will be the object of 
a future investigation. 

The purpose of this paper was to determine if quark 
cores with different properties would lead to different 
thermal evolution, in a way that one could use this 
information to constrain the microscopic properties 
of quark matter in hybrid stars. We have found a 
set of parameters that leads to a reasonably good 
agreement with the observed data of cooling neutron 
stars, although not all of them agree with the observed 
mass of J1614-2230. It is evident that there is room 
for improvement in the approach we used to model 
quark matter microscopically, and it is our intent to use 
the results and constraints put forth in this work, in 
addition to the latest observational data, to refine the 
current models for the EoS of quark matter. 



[1] F. Weber, Pulsars as astrophysical laboratories for nu- [2] N. K. Glendenning, Compact stars : nuclear physics, 
clear and particle physics, 1st ed. (Institute of Pliysics, particle physics, and general relativity, 1st ed. (Springer, 

Bristol, 1999) 2000) 



8 



F. Ozel, G. Baym, and T. Guver, arXiv/1002. 3153(2010) [29 
A. W. Steiner, J. Lattimer, and E. F. Brown, 
arXiv/1005.0811(2010) [30 
D. Page, J. Lattimer, M. Prakash, and A. W. Steiner, 
The Astrophysical Journal Supplement Series 155, 623 [31 
(2004) 

D. Page, J. Lattimer, M. Prakash, and A. W. Steiner, [32 
The Astrophysical Journal 707, 1131 (2009) 
D. Page, U. Geppert, and F. Weber, Nuclear Physics A [33 
777, 497 (2006) 

J. Lattimer, C. Pethick, M. Prakash, and P. Haensel, [34 
Physical review letters 66, 27012704 (1991) 
D. Blaschke, T. Klahn, and D. Voskresensky, The Astro- [35 
physical Journal 533, 406412 (2000) 
C Schaab, F. Weber, M. Weigel, & N. K. Glendenning, [36 
Nuclear Phys A 605, 531 (1996). [37 
K. P. Levenfish, & D. G. Yakovlev, Astronomy Letters 
20, 43 (1994). ' [38 

D. Page, M. Prakash, J. M. Lattimer, & A. W. Steiner, 
astro-ph/1110.5116 (2011). [39 

H. Grigorian, D. Blaschke, and D. Voskresensky, Physical 
Review C 71, 1 (2005) ' [40 

D. Blaschke, D. Voskresensky, and H. Grigorian, Nuclear 
Physics A 774, 815 (2006) [41 
P. Jaikumar, M. Prakash, and T. Schafer, Physical Re- [42 
view D 66, 1 (2002) 

I. Shovkovy and P. EUis, Physical Review C 66, 1 (2002) [43 
R. Tolman, Physical Review 55, 364 (1939) 
J. Oppenheimer and G. Volkoff, Physical Review 55, 374 [44 
(1939) 

H. Grigorian, D. Blaschke, and D. Aguilera, Physical Re- 
view C 69, 1 (2004) [45 
R. Negreiros, V. Dexheimer, and S. Schramm, Physical 
Review C 82, 1 (2010) [46 
V. A. Dexheimer and S. Schramm, Physical Review C 81 
(2010) [47 
N. K. Glendenning, Nuclear Physics A 493, 521 (1989) 
A. Chodos, R. Jaffe, K. Johnson, C. Thorn, and V. Weis- [ 
skopf. Physical Review D 9, 3471 (1974) 
A. Chodos, R. Jaffe, K. Johnson, and C. Thorn, Physical [49 
Review D 10, 2599 (1974) 

E. Farhi and R. Jaffe, Physical Review D 30, 2379 (1984) 

F. Weber, Progress in Particle and Nuclear Physics 54, [50 
193 (2005) 

C. Alcock, E. Farhi, and A. Olinto, The Astrophysical [51 

Journal 310, 261 (1986) 

V. Usov, Physical Review D 70, 14 (2004) 



R. Negreiros, F. Weber, M. Malheiro, and V. Usov, Phys- 
ical Review D 80, 083006 (2009) 

R. P. Negreiros, L Mishustin, S. Schramm, & F. Weber, 
Physical Review D 82, 103010 (2010). 
M. Alford, Annual Review of Nuclear and Particle Sci- 
ence 51, 131 (2001) 

M. Alford, A. Schmitt, K. Rajagopal, and T. Schafer, 
Reviews of Modern Physics 80, 1455 (2008) 
M. Alford, J. Bowers, and K. Rajagopal, Physical Review 
D 63, 074016 (2001) 

J. Bowers and K. Rajagopal, Physical Review D 66, 
065002 (2002) 

G. Lugones and J. E. Horvath, Physical Review D 66, 1 
(2002) 

M. Alford and S. Reddy, Physical Review D 67, 1 (2003) 
G. Baym, C. Pethick, and P. Sutherland, The Astrophys- 
ical Journal 170, 299 (1971) 

E. H. Gudmundsson, C. J. Pethick, and R. I. Epstein, 
The Astrophysical Journal 259, L19 (1982) 
R. L E. E. H. Gudmundsson, C. J. Pethick, The Astro- 
physical Journal 272, 286 (1983) 

D. Yakovlev, A. D. Kaminker, O. Y. Gnedin, and 
P. Haensel, Physics Reports 354, 1 (2001) 

N. Iwamoto, Annals of Physics 141, 1 (1982) 

E. Flowers and N. Itoh, The Astrophysical Journal 250, 
750 (1981) 

P. Haensel, Nuclear Physics B - Proceedings Supplements 
24, 23 (1991) 

C. M. Espinoza, a. G. Lyne, M. Kramer, R. N. Manch- 
ester, & V. M. Kaspi, The Astrophysical Journal 741, 
L13 (2011). 

A. Kurkela, P. Romatschke, A. Vuorinen, and B. Wu, 
arXiv/1006. 4062(2010) 

D. Blaschke, J. Berdermann, and R. Lastowiecki, 
arXiv/1009. 1181(2010) 

P. B. Demorest, T. Pennucci, S. M. Ransom, M. S. E. 
Roberts, & J. W. T. Hessels, Nature 467, 1081 (2010). 
D. Page, M. Prakash, J. Lattimer, & A. Steiner, Physical 
Review Letters 106, 081101 (2011). 

P. S. Shternin, D. G. Yakovlev, C. O. Heinke, W. C. G. 
Ho, & D. J. Patnaude, Monthly Notices of the Royal 
Astronomical Society; Letters 412, L108-L112 (2011). 
T. Noda, M. Hashimoto, Y. Matsuo et al, astro- 
ph/1109.1080 

R. Negreiros, S. Schreamm, & F. Weber, astro- 
ph/1103.3870 



